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Abstract — Constellation shaping is an energy-efficient strategy 
involving the transmission of lower-energy signals more fre- 
quently than higher-energy signals. Previous work has shown 
that shaping is particularly effective when used with coded 
amplitude phase-shift keying (APSK), a modulation that has been 
popularized recently due to its inclusion in the DVB-S2 standard. 
While shaped APSK can provide significant gains when used with 
standard off-the-shelf LDPC codes, such as the codes in the DVB- 
SI standard, additional non-negligible gains can be achieved by 
optimizing the LDPC code with respect to the shaped APSK 
modulation. In this paper, we optimize the degree distributions 
of the LDPC code used in conjunction with shaped APSK. The 
optimization process is an extension of the EXIT-chart technique 
of ten Brink, et al., which has been adapted to account for the 
shaped APSK modulation. We begin by constraining the code to 
have the same number of distinct variable-node degrees as the 
codes in the DVB-S2 standard, and show that the optimization 
provides 32- APSK systems with an additional coding gain of 0.34 
dB at a system rate of = 3 bits per symbol, compared to shaped 
systems that use the long LDPC code from the DVB-S2 standard. 
We then increase the number of allowed variable node degrees 
by one, and find that an additional 0.1 dB gain is achievable. 

L Introduction 

The combination of low-density parity-check (LDPC) codes 
and ampHtude phase- shift keying (APSK) has become popu- 
larized as a result of the widespread adoption of the Digital 
Video Broadcasting Satellite - Second Generation (DVB-S2) 
standard, which uses LDPC-coded APSK for the satellite 
delivery of high-definition television, interactive services, and 
data-content distribution |[T|. An APSK signal constellation 
consists of several concentric rings of symbols, with each 
ring having a different amplitude and containing symbols that 
are separated by a constant phase offset |2|. When combined 
with LDPC codes, APSK offers an attractive combination of 
spectral and energy efficiency, and is particularly well suited 
for the nonlinearities in satellite channels that arise due to the 
use of high-power amplifiers. 

In typical implementations, such as DVB-S2, the APSK 
symbols are selected with uniform probability. However, as 
we have shown in |[3|, |[4|, the performance of APSK may 
be significantly improved by selecting lower-energy signals 
more frequently than higher-energy signals. The strategy in 
|[3|, |[4j| is an instance of the general concept of constellation 
shaping |[5|, |[6|. Our specific strategy for shaping APSK 
is adapted from a strategy originally proposed by Le Goff, 
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et al., in |6| for bit-interleaved turbo-coded pulse-amplitude 
modulation (PAM). The strategy partitions the base signal 
constellation into two or more sub-constellations of increasing 
average energy and uses a shaping code |^ to select signals 
from the lower energy sub-constellations more often than 
the signals from the higher energy sub-constellations. While 
shaping increases the per-iteration complexity of the system, 
the error-rate performance can often be improved enough that 
fewer iterations are required compared to a uniform system 
operating at the same signal-to-noise ratio (SNR) |4^|. 

Shaping an APSK constellation can achieve significant 
gains, even if an "off-the-shelf" standard channel code is used. 
For instance, in | 3 1 we show a gain in excess of 0.6 dB when 
shaped APSK is combined with UMTS standardized turbo 
code at a rate of 3.85 bits/symbol, while in |4| shaping gains 
in excess of 0.7 dB was observed for LDPC-coded 3 2- APSK 
operating at a rate of 3 bits/symbol. While using standardized 
codes is convenient, a significant additional gain can be 
achieved by tailoring the code to match the characteristics of 
the shaped modulation. This fact can be visualized through 
the use of extrinsic information transfer (EXIT) charts ||7|. 
The EXIT chart of an LDPC-coded system comprises a pair 
of curves, which need be properly matched for optimal perfor- 
mance. Since shaping changes the corresponding EXIT curve, 
a code whose EXIT curves are well matched in a uniform 
modulation scheme may not yield the optimal performance 
for a shaped modulation system. By changing the variable- 
node degree distribution, it is possible to identify a code for 
the shaped system with better matched EXIT curves. As we 
show, the improvement in performance can be significant. 

The remainder of this paper is organized as follows. Section 
[Il| provides a model for bit-interleaved LDPC-coded APSK 
with constellation shaping. Section III compares the perfor- 
mance of a shaped system against that of a comparable 
uniform system operating at the same rate, with both systems 



using codes from the DVB-S2 standard. Section |IV] uses EXIT 
charts (Tj to identify better degree distributions for the uniform 
APSK system, while section [V| extends the technique to find 
optimal degree distributions for the shaped system. In these 
sections, the optimal degree distributions are given along with 
simulation results for LDPC-coded 32-APSK operating over 
an additive white Gaussian noise (AWGN) channel at a rate 
of 3 bits/symbol with LDPC codes of two lengths. Finally, the 



paper concludes in Section VI 
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Fig. 1. Transmitter structure. 



II. System Model 



The system's transmitter is shown in Fig. [T] A length-/cc 
vector b of information bits is fed into a rate Rc = kc/ric 
LDPC encoder, which produces a length-nc codeword u. 
The first interleaver Hi produces a vector v after random 
permutation. A bit separator separates v into m streams, 
where m = log2 M and M is the number of symbols in 
the APSK constellation. For the remainder of this paper, we 
assume that M — 32, and hence, m = 5. The first stream, 
denoted d, is passed through a rate Rs = kg/ rig shaping 
encoder to produce the shaped codeword c, which is then 
permuted by a second interleaver 112 to produce the vector si. 
The purpose of the shaping encoder is to produce a nonuniform 
output with a probability po of that is higher than the 
probability pi of 1. The detailed operation of the shaping 
encoder is described in |4|. The remaining streams at the 
output of the bit separator are denoted {s2, S5} and are of 
the same length as Si. The streams flow into a bit combiner, 
which passes groups of m = 5 bits to the APSK modulator. 

The 3 2- APSK constellation is shown in Fig. [2] and follows 
the DVB-S2 standard, except the two most significant bits 
(MSBs) are exchanged to allow the shaping bit to be the MSB. 
The bit from Si is called the shaping bit and it partitions the 
constellation into two sub-constellations. The first partition, 
which is labeled with an MSB equal to 0, is selected with 
probability pq, and contains the signals in the inner two rings. 
The second partition, which is labeled with an MSB equal 
to 1, is selected with probability pi, and contains the signals 
in just the outer ring. Because the shaping code ensures that 
Po > Pi, the signals in the inner two rings are more likely to 
be selected than the high-energy signals of the outer ring. 

The overall system rate R is the number of information 
bits per modulated symbol, and is related to the rate Rc of 
the LDPC code and the rate Rg of the shaping code by R = 
Rc{m-\- Rs — 1). When uniform modulation is used, Rs = I 
and there is no shaping code. It is important to note that the 
choice of Rc and Rs was made to provide a fair comparison 
between the uniform and shaped system at the same R, so the 
shaped system must use an LDPC code with a higher rate Rc 
than is used by the uniform system. If the LDPC codes are 
selected from among those in the DVB-S2 standard (which has 
just eleven different rates), then the choice of Rs is limited 
and may be suboptimal. Indeed, one of our motivations for 
designing new LDPC codes for the shaped system was so that 
we could select an LDPC code rate Rc that can be matched 
with a shaping code with a better rate Rs. 
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Fig. 2. 32- APSK constellation. The shaping bit is the MSB. 

III. Baseline Unoptimized Codes 

Before attempting to optimize the LDPC code, we compare 
the performance of a uniform LDPC-coded system against 
a comparable shaped system, but restrict the LDPC codes 
to come from the DVB-S2 standard (H in both cases. This 
will serve as a baseline, against which we can compare our 
optimized codes. For both shaped and uniform systems, the 
overall rate is = 3 bits/symbol. To achieve this rate with 32- 
APSK, the uniform system must use a rate Rc — 3/5 LDPC 
code, while for the shaped system, a = 2/3 LDPC code 
combined with a = 1/2 shaping code is used for the overall 
rate to remain R = 3 bits/symbol. The shaping code with 
{usjks) = (4,2) provides a reasonably good po while still 
allowing rates Rc from the DVB-S2 standard to be used (in 
Section |Vl we will consider custom-rate LDPC codes that will 
allow a shaping code with better po to be used). The codewords 
of this shaping code are: {(0000), (0100), (0010), (0001)}. 

The uniform and shaped systems were simulated over an 
AWGN channel using long (length ric = 64,800) LDPC 
codes. The APSK modulation was configured using ring- 
radii ratios that maximize capacit}Q The receiver for the 
uniform system uses the bit-interleaved modulation with itera- 
tive decoding (BICM-ID) architecture described in jSj, which 
involves the iterative exchange of information between the 

^Let 71 be the ratio of the radius of the middle ring of the 32- APSK 
constellation to the radius of the inner ring, and 72 be the ratio of the radius 
of the outer ring to the radius of the inner ring. Together, the couple {71 , 72} 
constitutes the set of ring-radii ratios. For both shaped and uniform 32-APSK, 
the ring-radii ratios that maximize capacity are (2.64,4.64) |3j 
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Uniform with optimized 
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Uniform with optimized 
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Fig. 3. Bit-error rate of 32-APSK in AWGN at rate R = 3 bits/symbol with a 
length Tic = 64 800 LDPC code. The soHd lines are for the uniform system; 
the rightmost curve shows the performance of uniform modulation with a 
standard DVB-S2 code, while the two curve to its left show the performance 
of uniform modulation with optimized D = 3 and = 4 LDPC codes, 
respectively. The dashed lines are for the shaped system; with the left two 
curves show the performance of shaped modulation with optimized D = 3 
and D = 4 LDPC codes, and the third curve from left shows the performance 
of shaped modulation with a standard DVB-S2 code. 

APSK demapper and the LDPC decoder. A similar BICM- 
ID receiver was adopted for the shaped system, but with the 
shaping decoder included in the overall feedback loop. See 
(41 for details of the receiver operation. For both systems, a 
maximum of 100 decoder iterations were performed, though 
the decoder was halted if all the parity-checks were satisfied. 

In addition to other curves that will be described later, Fig.|3] 
shows the simulated bit-error rates of the two cases described 
in this section: (1) the uniform system with the rate Rc = 
3/5 DVB-S2 standardized LDPC code and BICM-ID reception 
(rightmost curve); and (2) the shaped system using the Rc = 
2/3 DVB-S2 standardized LDPC code, Rs = 2/4 shaping 
code, and iterative receiver of |4| (third curve from left). The 
LDPC codes that are marked "optimized" will be described 
in sections |IV] and |V] From Fig. [3] it is observed that the 
uniform and shaped systems require E^/Nq equal to 5.42 dB 
and 4.96 dB to achieve a BER equal to 10~^, respectively. 
Thus, the shaped system achieves a gain of 0.46 dB relative 
to the uniform system when both use BICM-ID reception. The 
capacity limit of 32-APSK with po = 0.8125, corresponding 
to Rs = 2/4, is at ^/TVq = 3.83 dB ji), thus this shaped 
system is 1.13 dB away from the capacity. 

IV. Code Optimization for Uniform Modulation 

We now optimize the LDPC code. Our goal is to adapt 
the EXIT-chart technique of |7 1 to account for the behavior of 
the shaped modulation. However, before optimizing the LDPC 
code used by the shaped system, we begin by reviewing how 
the EXIT charts can be used to optimize the LDPC code in 
the system with uniform modulation. To facilitate systematic 
encoding, we limit the LDPC code to be an extended irregular 
repeat accumulate (eIRA) code |9 |, which is also the code 



structure defined in DVB-S2 standard. The parity-check matrix 
H of such a code contains an arbitrary matrix for the first kc 
columns and a dual-diagonal matrix as its last ric — kc columns. 
The code is further limited to be check regular, and its check- 
node degree dc is kept the same as the value of the DVB-S2 
standardized code of the same rate and length. For instance, 
the standardized rate Rc = 3/5 LDPC code corresponds to 
the parity-check matrix with dimension (25 920 x 64 800) and 
constant row weight dc = 11, and the comparable optimized 
code is constrained to have the same check-node degree. 

Since the check-node degree is fixed, the optimization prob- 
lem is to select the optimal variable-node degree distribution. 
Following 1 7 1, the variable-node degree is indicated by dy^i, 
{i = 1,...,!)). The fraction of nodes having degree dy^i is 
denoted by a^, while the fraction of edges incident to variable 
nodes of degree dy^i is denoted by bi. First we limit the number 
of distinct variable-node degrees to be I) = 3, since this is 
the same number allowed in the DVB-S2 standard. We then 
increase D by one, and optimize LDPC codes with = 4 
using the same approach. Due to the dual-diagonal part of the 
H matrix, the minimunj^ variable-node degree is dy^i = 2. 
This implies that the fraction of degree-2 nodes is (at leasj^ 
ai = (n — k)/n; i.e., the fraction of columns corresponding 
to the dual-diagonal part of H. To constrain complexity, the 
maximum node degree is limited to be no greater than 25; 
i.e., dy^D < 25. When D = 3 is considered, the key step 
in the optimization is to determine the integer values of dy^2 
and dy^s. These are related to the constraints on the code rate 
and total number of Is (set by the check-node degree) in the 
parity-check matrix, thus dy^s can be determined once dy^2 
is set. Furthermore, once dy^2 and dy^s are determined, the 
corresponding degree distributions, a2 and as, can be found by 
using the check-node degree dc. The optimization is therefore 
over just one parameter. On the other hand, the optimization is 
over two parameters when = 4; two of the values in the 3- 
tuple {dy^2^dy^s^dy^4} must first be fixed before the remaining 
value can be determined. 

An EXIT chart comprises two EXIT curves |[7|. One EXIT 
curve corresponds to the variable-node decoder (VND), and 
the other one corresponds to the check-node decoder (CND). 
The EXIT-chart based LDPC code optimization fits the two 
EXIT curves together by picking appropriate variable-node 
degrees. The VND curve characterizes both the variable nodes 
of the LDPC code and the modulation part. It is created by first 
generating a transfer characteristic for just the modulation and 
its detector at the given Ss/Nq. This detector characteristic 
is denoted by /£;,det(^a, ^s/^o), where I a is computed 
between the demodulator's a priori inputs and the modulator's 
input bits. For APSK, /£;,det cannot be generated in closed 
form for nonzero I a, and thus it is generated through Monte 
Carlo simulation under the assumption that the demodulator's 

^In DVB-S2, there is a single degree- 1 variable node corresponding to the 
last column of the H matrix. While our code design maintains this degree- 1 
node, we do not explicitly list the fraction of nodes connected to degree- 1 
nodes since it will approach zero for large ric. 

^We found no advantage to using ai > (n — k)/n. 



a priori input is conditionally Gaussian. 

Having found the detector characteristic, the VND curve for 
degree-d^ nodes is found by using 

^£;,VND {lA^dy^Ss/No) = 

J (^^K - W-HIa)^] + [J-\Ie,bet{Ia.Ss/NoW^ 

(1) 

where the J-function is given in {T\ and can be computed 
using the truncated series representation of 1 10 |. Note that ([T]) 
only gives the EXIT curve for a single variable-node degree 
dy, and therefore represents the VND curve for a regular code. 
The VND curve for an irregular LDPC code is expressed as: 

D 

,VNB {lA,dy^i,Es/No) . (2) 

The CND curve is given in |7 | as: 

/£;,cnd(/a,4) = l-j(^Vdc-l-J-\l-lA))0) 

where I a is the mutual information at the input of the check 
nodes. The EXIT chart is drawn by noting that the /£;,cnd 
produced by the CND becomes the Ia at the input to the 
VND (which we denote /a,vnd), while the /£;,vnd produced 
by the VND becomes the I a at the input to the CND (which 
we denote /a,cnd)- The chart plots the VND and CND curves 
with /a, VND = ^£;,CND on the horizontal axis and /£;,vnd = 
^A,CND on the vertical axis. For a given degree distribution, 
the threshold is the value of Ss/Nq for which the VND and 
CND curve just barely touch. 

For I) = 3, all feasible combinations of {dy^2jdy^s} were 
considered and those that lead to low EXIT-chart thresholds 
were identified. One code found to be better than rate Rc = 
3/5 DVB-S2 standard code when used with 32-APSK has the 
following degree distribution: 

dy^i =2 ai= 0.40 bi 0.182 
dy[2 =4 a2 = 0.52 62 0.473 
dy^s = 19 ^3 = 0.08 bs = 0.345 

The process was repeated for D = A. Because the optimiza- 
tion for D = 4 was over two parameters, there were several 
degree distributions offering low EXIT-chart thresholds. A 
total of 50 candidates were identified, and the realized code 
for each was simulated. The code offering the best simulated 
performance has the following degree distribution: 
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The uniform 32-APSK system was simulated using a length 
nc = 64, 800 code with these two degree distributions, and the 
resulting BER curves with BICM-ID reception are shown in 
Fig. |3] The second and third curves from right are generated 
by using optimized rate Rc = 3/5 LDPC code with D = 3 
and D = 4. These codes indicate the gain of 0.14 dB and 0.30 
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Fig. 4. EXIT chart for shaped 32-APSK operating at system rate R = 3 in 
AWGN at Sb/No = 4.73 dB. Two different VND curves are shown, one that 
is drawn from the DVB-S2 standard and the other is optimized with D=4. 

dB relative to the DVB-S2 standard code at a BER equal to 
10~^. The LDPC code with one more variable-node degree 
(D = 4) achieves better performance. 

V. Code Optimization for Shaped Modulation 

When shaping is used, its effect is absorbed into the 
VND curve. Whereas generating the detector characteristic for 
the uniform case involves measuring the mutual information 
between the modulator input and demodulator output, the 
detector characteristic for the non-uniform constellation must 
account for the shaping code. In reference to Fig. [T] a bit 
sequence v is randomly generated and passed through the 
pictured processing to produce the modulated sequence x. The 
received noisy symbols are processed to produce the extrinsic 
information Le(v), and the detector characteristic is found by 
computing the mutual information between v and Le(v). As 
an example. Fig. [4] shows the VND curve of the DVB-S2 
standard rate Rc = 2/3 LDPC code combined with a (4,2) 
shaping code. Also shown is the CND curve for check-node 
degree dc = 10, which is the same dc of the rate Rc = 2/3 
LDPC code in the DVB-S2 standard. 

In the previous examples, both the LDPC code rate Rc and 
the overall rate R are chosen to be compatible with the code 
rates defined in the DVB-S2 standard. However, due to the 
limited number of LDPC code rates, the corresponding set 
of possible shaping code rates may lead to suboptimal po. 
For instance, the overall system rate R = 3 bits/symbol can 
be achieved with an LDPC code rate of Rc = 2/3 and a 
shaping code rate of Rs = 1/2. Using the (us^ks) = (4,2) 
shaping code will satisfy this requirement, but the resulting 
value of po is 0.8125, which differs from the optimal value 
of po = 0.716 identified in js]. A (3, 2) shaping code, with 
codewords {(000), (100), (010), (001)}, gives a po = 0.75, 
which is closer to the optimal, but requires LDPC code rate 
Rc = 9/14 to achieve an overall rate of R = 3 bits/symbol. 

There is no rate Rc = 9/14 code defined in DVB-S2 
standard. However, a new code can be designed for this rate 
by using EXIT-chart based techniques. As with the rate-2/3 
standardized codes, the check node degree of rate Rc = 9/14 
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Fig. 5. Bit-error rate of 32-APSK in AWGN at rate R = 3 bits/symbol with 
short length LDPC codes (ric = 16 200). The soHd lines are for the uniform 
system; from right to left, the curves show uniform system using Rc = 3/5 
DVB-S2 standard code, uniform system with optimized code with D = 3 
and D = 4. The dashed lines are for the shaped system; from left to right, 
the curves show performance using optimized Rc = 9/14 LDPC code with 
D = 4 and D = 3 combined with (3, 2) shaping code, the third curve from 
left shows the performance of standard Rc = 2/3 code with (4, 2) shaping 
code 

code is set to be = 10. For D = 3, the code found to have 
the best performance has variable-node degree distribution: 

dy^i =2 ai= 0.357 bi = 0.200 
dy[2 = 3 a2 = 0.558 62 = 0.469 
dy^s = 14 as= 0.085 63 = 0.331 

For I) = 4, the code found to have the best performance 
has variable-node degree distribution: 
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The VND curve for this code along with the (3, 2) shaping 
code is shown in Fig. |4] The difference compared with the 
standard LDPC code is subtle, but we can still observe an 
open tunnel for the optimized LDPC code. 

The BER performance of the two Rc = 9/14 length Uc = 
64 800 LDPC codes combined with a i^^ = 2/3 shaping code 
are shown in Fig.|3] The leftmost curve is generated by D = 4 
LDPC code, and the curve to its right is generated by using 
D = 3 LDPC code. The results indicate a gain at BER 10"^ of 
0.44 dB compared to the shaped system that uses the DVB- 
S2 standard rate Rc = 2/3 LDPC code along with a (4,2) 
shaping code, about 0.1 dB of the gain can be attributed to 
better designed LDPC code with more variable node degree. 

While the results in Fig. |3] are for long (n^ = 64 800) LDPC 
codes. Fig. |5] compares the bit-error performance of shorter 
length nc = 16 200 LDPC codes. In particular, the three solid 
curves on the right show the performance of uniform systems 
using short Rc = 3/b LDPC codes, with the second and third 
curve from right corresponding to the optimized LDPC code 
with D = 3 and D = A, respectively. The three dashed curves 
on the left correspond to shaped systems. The third curve from 



left is generated using the = 2/3 LDPC code from DVB- 
S2 standard combined with a (us^ks) = (4,2) shaping code. 
The first and second curve from left incorporate the short Rc = 
9/14 LDPC code optimized for shaping system with = 4 
and D = 3, respectively. At a BER of 10~^, optimizing the 
code for the uniform system provides a gain over the DVB-S2 
standardized code of 0.37 dB, while optimizing the code for 
the shaped system provides a gain of 0.43 dB. 

VL Conclusion 

Constellation shaping is particularly effective for APSK 
modulation. Even if "off-the-shelf" codes, such as those in the 
DVB-S2 standard, are used with shaping, impressive gains can 
be achieved. However, an additional and significant gain can 
be achieved by optimizing the code with respect to the shaped 
modulation. In the case of LDPC-coded 32-APSK operating at 
rate R = 3, optimizing the degree distribution while allowing 
D = 4 distinct variable degrees provides an additional gain 
of 0.44 dB for long (n^ = 64, 800) LDPC codes and 0.43 
dB for short (Uc = 16, 200) LDPC codes compared to shaped 
systems that are constrained to use codes from the DVB-S2 
standard. If a shaped 32-APSK system with optimized LDPC 
codes is compared against a uniform system of the same rate 
(R = 3), the shaping gain is between 0.6 dB and 0.73 dB, 
depending on the length of the code and number of distinct 
variable node degrees (D). Altogether, the combination of 
shaping and LDPC degree optimization has the potential to 
provide a combined (shaping+coding) gain in excess of 1 dB 
compared to the DVB-S2 standard, which neither shapes the 
modulation nor optimizes the LDPC degree distribution, and 
within 0.7 dB of the capacity limit for 32-APSK. 
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